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ABSTRACT 

The Hubble radius is a particular manifestation of the Universe's gravitational hori- 
zon, Rh(h)) = c/Hq, the distance beyond which physical processes remain unob- 
servable to us at the present epoch. Based on recent observations of the cosmic mi- 
crowave background (CMB) with WMAP, and ground-based and HST searches for 
Type la supernovae, we now know that Rh(to) ~ 1 3.5 Glyr. This coincides with the 
maximum distance (cto ~ 13.7 Glyr) light could have traveled since the big bang. 
However, the physical meaning of is still not universally understood or accepted, 
though the minimalist view holds that it is merely the proper distance at which the 
rate of cosmic recession reaches the speed of light c. Even so, it is sometimes argued 
that we can see light from sources beyond R^, the claim being that R^ lies at a redshift 
of only ~2, whereas the CMB was produced at a much greater redshift (~1 100). In 
this paper, we build on recent developments with the gravitational radius by actually 
calculating null geodesies for a broad range of FRW cosmologies, to show — at least 
in the specific cases we consider here, including ACDM — that no photon trajectories 
reaching us today could have ever crossed Rh(to)- We therefore confirm that the cur- 
rent Hubble radius, contrary to a commonly held misconception, is indeed the limit 
to our observability. We find that the size of the visible universe in ACDM, measured 
as a proper distance, is approximately 0.45cfo. 
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1 INTRODUCTION 

The standard model of cosmology is confronted with several unpalatable coincidences, suggesting 
that we do not yet have a fully consistent picture of the Universe's dynamical expansion (see, 
e.g., Melia & Shevchuk 2011). Part of the problem is that cosmological observations can only be 
interpreted from within the context of a pre-assumed model, and the data can be quite compliant. 

The ACDM (Cold Dark Matter with a cosmological constant A) model has been without peer 
in cosmology (see, e.g., Spergel et al. 2003, and Tegmark et al. 2004). For example, this model 
has been used with complementary measurements of the cosmic microwave background (CMB) 
radiation to determine that the Universe is flat, so its energy density p is at (or very near) the 
"critical" density p c = 3c 2 H 2 /SnG. But among the many peculiarities of this description of the 
universe is the inference, based on current observations, that the density p d of dark energy must 
itself be of order p c . (In the context of ACDM, the best fit to the WMAP data indicates that dark 
energy represents approximately 73% of the total p « p c ; see Spergel et al. 2003.) Dark energy is 
often thought to be the manifestation of the aforementioned cosmological constant, A, though no 
reasonable explanation has yet been offered as to why such a fixed, universal density ought to exist 
at this scale. It is well known that if A is associated with the energy of the vacuum in quantum 
theory, it should have a scale representative of phase transitions in the early Universe — 120 orders 
of magnitude greater than p c . 

Many workers have attempted to circumvent these difficulties by proposing alternative forms 
of dark energy, including Quintessence (Ratra & Peebles 1988; Wetterich 1988), which represents 
an evolving canonical scalar field with an inflation-inducing potential, a Chameleon field (see, e.g., 
Mota & Barrow 2004; Khoury & Weltman 2004; Brax et al. 2004) in which the scalar field couples 
to the baryon energy density and varies from solar system to cosmological scales, and modified 
gravity, arising out of both string motivated, or General Relativity modified actions (Capozziello et 
al. 2003; Nojiri & Odintsov 2003; Carroll et al. 2004), which introduce large length scale correc- 
tions modifying the late time evolution of the Universe. The actual number of suggested remedies 
is far greater than this small, illustrative sample. 

An equally perplexing puzzle with ACDM has been dubbed the "coincidence problem," arising 
from the peculiar near-simultaneous convergence of the matter energy density p m and pd towards 
p c in the present epoch. Though p m and p d are expected to change at different rates as the Universe 
expands (particularly if dark energy is a cosmological constant) they are nearly equal in the present 
epoch, implying that we live at a special time in cosmic history. 
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In a recent paper (Melia & Shevchuk 201 1), we proposed an explanation for yet another dis- 
turbing coincidence, having to do with the apparent equality of our gravitational (or Hubble) radius 
R h with the distance cto light could have traveled since the big bang (in terms of the presumed cur- 
rent age t of the Universe). This equality has received some scrutiny in recent years (Melia 2003, 
2007, 2009, Melia & Abdelqader 2009, van Oirschot et al. 2010; see also Lima 2007 for a related, 
though unpublished, work). 

Unfortunately, there is still some confusion regarding the properties of Rh due to a misun- 
derstanding of the role it plays in our observations. For example, it is sometimes suggested (see, 
e.g., Davis & Lineweaver 2004; van Oirschot et al. 2010) that sources beyond Rh(k) are observ- 
able today, which is certainly not the case. We will therefore begin by elaborating upon what the 
gravitational radius R h is — and what it is not. Though first defined in Melia (2007), an unrecog- 
nized form of Rh actually appeared in de Sitter's (1917) own account of his spacetime metric. 
And we will advance the discussion further by actually calculating photon trajectories for vari- 
ous well-studied Friedmann-Robertson- Walker (FRW) Cosmologies, demonstrating that the null 
geodesies reaching us at t have never crossed Rh(t ). Some come close, and in one case — the de 
Sitter model — they approach R h asymptotically as t recedes to our infinite past. Our conclusion in 
this paper will be that Rh(to) is a real limit to our observability at the present time t . 



2 THE GRAVITATIONAL (HUBBLE) RADIUS 

Standard cosmology is based on the Robertson- Walker metric for a spatially homogeneous and 
isotropic three-dimensional space. In terms of the proper time t measured by a comoving observer, 
and the corresponding radial (r) and angular coordinates (6 and <p) in the comoving frame, an 
interval ds in this metric is written as 

ds 2 = c 2 dt 2 - a 2 (t)[dr 2 {\ - kr 2 y l + r 2 {d6 2 + sin 2 6d(p 2 )] , (1) 

where a{t) is the expansion factor and the constant A: is +1 for a closed universe, for a flat, open 
universe, or - 1 for an open universe. 

In recent work (Melia 2007, 2009; Melia & Abdelqader 2009), we demonstrated the useful- 
ness of examining properties of the metric in terms of both co-moving coordinates {ct, r, 6, (p) and 
observer-dependent coordinates {cT,R, 9, 0), where R is the so-called proper radius a{t)r. Whereas 
{ct, r, 9, (p) describe events in a frame "falling" freely with the cosmic expansion, the second set 
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of coordinates are referenced to a particular individual who describes the spacetime relative to the 
origin at his location 3 

For a flat universe (k = 0), it is straightforward to show that with this definition of R, Equa- 
tion (1) becomes 

r / p \ I 2 

-Q^dR 2 -R 2 dQ 2 , (2) 



ds 2 = O 



Rv 



where the function 

,2 



l*h, 

signals the dependence of the metric on the proximity of the observation radius R to the gravita- 
tional radius R h . As shown in Melia & Abdelqader (2009), the exact form of Rh depends on the 
constituents of the Universe. For example, in de Sitter space (which contains only a cosmological 
constant A), R h = c/H , in terms of the (time-invariant) Hubble constant H . It is equal to let in a 
radiation-dominated universe and 3ct/2 when the Universe contains only (visible and dark) matter. 

Quite generally, R h is the radius at which a proper sphere encloses sufficient mass-energy to 
turn it into a Schwarzschild surface for an observer at the origin of the coordinates (see Melia & 
Shevchuk 201 1). That is, Rh is given in all cases by the expression 

(4) 

where M(R h ) = (4tt/3)^p/c 2 , in terms of the total energy density p. Thus, R h = (3c 4 /8nGp) l/2 
which, for a flat universe, may also be written more simply as Rh(t) = c/H Q (t). 

Although not defined in this fashion, the Hubble radius c/H (t) (more commonly encountered 
when co-moving coordinates are used) is therefore seen to coincide with the gravitational radius 
R h emerging directly from the Robertson- Walker metric written in terms of R. What this means, 
of course, is that the Hubble radius is not a mere empirical artifact of the expanding universe, but 
actually represents the radius at which a sphere centered on the observer contains sufficient mass- 
energy for its surface to function as a static horizon. Of course, that also means that the speed of 
expansion a distance R h away from us must be equal to c, just as the speed of matter falling into a 
compact object reaches c at the black hole's horizon (and this is in fact the criterion used to define 
the Hubble radius in the first place). 

The reason why R^ is an essential ingredient of the metric written in the form of Equation (2) 
can be understood in the context of the Birkhoff (1923) theorem and its corollary (see Melia 2007). 



1 It is worth mentioning that prior to the introduction of co-moving coordinates in the 1920's, the cosmological spacetime was actually conven- 
tionally expressed in terms of these observer-dependent coordinates (see, e.g., de Sitter 1917.) 
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The latter states that the metric inside an empty spherical cavity, at the center of a spherically 
symmetric system, must be equivalent to the flat-space Minkowski metric. Space must be flat in a 
spherical cavity even if the system is infinite. Thus, the metric between the edge of a cavity and 
a spherically symmetric mass placed at its center is necessarily of the Schwarzschild type. The 
worldlines linked to an observer in this region are curved relative to the center of the cavity in 
a manner determined solely by the mass we have placed there. The implication is that any two 
points within a medium with non-zero p experience a net acceleration (or deceleration) toward (or 
away) from each other, determined solely by how much mass-energy is contained between them. 
This is why, of course, the Universe cannot be static, a concept that eluded Einstein himself, since 
his thinking on this subject preceded that of Birkhoff. 

Yet the physical meaning of R h is still elusive to many, possibly because of the widely held 
belief that all horizons must necessarily be asymptotic surfaces attained when t —> <x>. The so- 
called event horizon (see Rindler 1956) is indeed of this type, representing the ultimate limit of 
our observability to the end of time. However, R h is not in this category — nor should it be. Unlike 
the event horizon, the gravitational radius is a time-dependent quantity that increases in value 
at a rate determined by the evolving constituents of the Universe, specifically, the value of the 
equation-of-state parameter w, defined by the relation p = wp, where p is the pressure and p is the 
energy density. For some cosmologies, Rh may turn into the event horizon when the cosmic time 
approaches infinity. 

Note from the definition of R h that Rh/Rh = ~p/2p, so with Equation (4) 

R h = 3 -{\+w)c. (5) 

Clearly, R h is constant only for the de Sitter metric, where w = -1 and therefore Rh = 0. For 
all other values of w > -1, Rh > 0. So in ACDM, for example, where the Universe is currently 
dominated by a blend of matter and dark energy, Rh > 0. If dark energy were a cosmological 
constant, however, the Universe would eventually become de Sitter as the density of matter drops 
to zero, and we would therefore expect R h to then asymptotically approach a constant value equal 
to the radius of the event horizon in ACDM. 

What this means, then, is that the current location of R h affects what we can observe right 
now, at time t since the big bang; it is not — and is not meant to be — an indication of how far 
we will see in our future. One must therefore be careful interpreting spacetime diagrams, such 
as Figure 3 in van Oirschot et al. (2010). Consider that for an Einstein-de Sitter universe (with 
w = 0, the value used to construct this figure), Rh = (3/2)c so our observations never encounter 
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R h because no matter how far light travels, it never overtakes R h . But that doesn't mean that this 
must always be the case. In fact, Rh < c for any cosmology with w < -1/3 — which is presumably 
the situation with ACDM. One of the principal goals of this paper is to calculate the actual photon 
trajectories for well-known FRW cosmologies, such as ACDM, in order to demonstrate how and 
why R h functions as a limit to our observability today (at time t ). 



3 NULL GEODESICS IN FRW COSMOLOGIES 

Let us now derive the equation describing photon trajectories in a cosmology consistent with the 
FRW metric (Equation 1). From the definition of proper radius, we see that 

R = dr + ar . (6) 

But the null condition for geodesies (see, e.g., Ellis & Rothman 1993) applied to Equation (1) 
yields 

dr 

cdt = -a(t)—= , (7) 
Vl - kr 2 

where we have assumed propagation of the photon towards the origin along a radius. The best 
indications we have today are that the universe is flat so, for simplicity, we will assume k = in 
all the calculations described below, and therefore (for a photon approaching us) r = -c/a. Thus, 
we can write Equation (6) as follows: 

*r-(|-l). W 
where we have added a subscript y to emphasize the fact that this represents the proper radius of 
a photon propagating towards us. Note that in this expression, both R y and R h are functions of 
cosmic time t. The gravitational radius must therefore be calculated according to Equation (5). 

The analysis of the WMAP data, within the context of the standard model, points to an age for 
the Universe of to ~ 13.7 billion years. But we also know p quite accurately now, and we can use it 
to calculate Rh(t ), which appears to be approximately 13.5 billion light-years, equal to ct within 
the measurement errors. This is the unlikely coincidence we alluded to in the introduction, because 
there clearly is no particular reason why R h in Equation (5) should be equal to c, especially if w 
changes with time. 

The implausibility of this equality and its possible interpretation have been discussed elsewhere 
(e.g., Melia & Shevchuk 2011), so we will not dwell on them here. However, we will use the 
inferred value of R h today as one of our boundary conditions. In principle, Equation (8) may be 
integrated either forwards or backwards, but in reality, we are more familiar with the physical 
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conditions now, at time t , so it makes practical sense to begin the solution of this equation at 
t = to, where R y = and Rh(to) = ct . 

Not surprisingly, R y equals -c today, since this must represent the actual speed of light mea- 
sured in our local frame at the origin of our coordinates. Notice, however, that R y differs from -c 
away from the observer. This type of phenomenon often gives rise to misinterpretation and con- 
fusion, stemming from the fact that many consider R to be the actual physical distance measured 
by a single observer, and that dR/dt is therefore the physical speed. But in general relativity, the 
proper velocity measured by an individual is actually v = -\Jg RR /g TT dR/dT, calculated in terms of 
the proper distance R and proper time T in his/her frame, where the metric coefficients g RR and g TT 
are independent of T. Light satisfies the null condition ds = 0, and therefore v is always equal to 
c, regardless of whether the frame of reference is inertial or not. Written in terms of the co-moving 
coordinates, the metric in Equation (1) does not satisfy these conditions, and so R y is generally not 
equal to c. 

This happens because t only represents the physical time on the clock at the observer's location 
(in other words, t = T only at R = 0). If this observer were to measure T at any other position, 
(s)he would infer a value different from the local proper time t at that radius. And because t is used 
as a common time everywhere, the quantity R(t) at best represents the sum of all the incremental 
segments measured by a collaboration of observers — each at time t in their own rest frame — 
stretched out from one spacetime point to another (see, e.g., Weinberg 1972). Thus, although R(t) 
is the "proper" distance defined in terms of the cosmic time t, it actually does not represent the 
physical distance between these two points for an individual observer whose time coordinate is 
T in his/her frame. A better description for R would be that it represents a "community" distance 
between two spacetime points and, as a result, R y is not constrained to have the value -c away 
from the origin. 

In the next section, we will discuss the solutions to Equation (8) for various assumed cosmolo- 
gies, including ACDM, the current standard model. 

4 RESULTS AND DISCUSSION 

The de Sitter universe has no ordinary matter or radiation. Its dynamics is dictated solely by a 
cosmological constant A, which results in an expansion factor a(t) = e Hat . As we noted earlier 
(following Equation 3), the gravitational horizon for this spacetime is constant R h = c/H and 
w = -1. By design, the energy density p is also constant. A plot of R y , in units of R h (t ), is shown 
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Figure 2. Same as figure 1, except here for the Einstein-de Sitter universe, consisting entirely of matter. 

as a function of t/t in figure 1. In de Sitter, the maximum proper distance of photons reaching us 
today was « 0.65i? n (/o), which occurred at cosmic time t = 0. Never along its geodesic path did 
light's proper distance from us exceed our gravitational horizon today. 

This feature is echoed by the Einstein-de Sitter universe, which consists entirely of matter. The 
cosmological constant and radiation are both absent, so w = 0. In this case, the scale factor grows 
as a(t) = (3H t/2) 2/3 (see, e.g., Melia & Abdelqader 2009), so the density is infinite at t = and 
decreases monotonically as t —> t . The gravitational radius (R h = 3c//2) grows linearly with time, 
but is always greater than ct. The null geodesic for this cosmology is shown in figure 2, where we 
see that R y attains its maximum value ~0.3/? n (fo) about 2.6 Gyr after the big bang. 

When the universal expansion is driven solely by radiation, w = +1/3 and R h = let. The 
expansion factor is then given by a(t) = (2H Q t) 1/2 , and the corresponding null geodesic is shown 
in figure 3. The proper distance R y of light attains its maximum value of ~0.4/? h (? ) at t « 0At . 

For a universe dominated by a combination of matter and radiation, the expansion factor also 
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Figure 3. Same as figure 1, except for a universe dominated by radiation, for which w = +1/3 and Sj, = let. 
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Figure 4. Same as figure 1 , except for a universe dominated by matter and radiation. 



grows monotonically with cosmic time t, though the value of w falls sharply from about +1 /3 in the 
early universe (where radiation dominates) and tapers off towards zero as we approach the present 
time, where the effect of radiation is relatively insignificant. The photon geodesic for this case is 
shown in figure 4, where we see it attaining its maximum proper distance from us, ~0.37? h (? ), 
roughly 3 Gyr after the big bang. As with the other cases we've considered, the photon's trajectory 
never takes it beyond our current gravitational horizon R h (to). 

In ACDM, the standard model of cosmology, the density p is comprised of three principal 
components, 

P = Pm + Pr + Pde (9) 

where, following convention, we designate the matter, radiation, and dark energy densities, re- 
spectively, as p m , p r , and p de . We will also assume that these energy densities scale according to 
p m oc 6f 3 , p r oc a' 4 , and p de oc f(a). (If dark energy is indeed a cosmological constant A, which we 
assume for this calculation, then f(a) = constant.) In this case, there is no simple analytical solu- 
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0.0 0.2 0.4 0.6 0.8 1.0 

t/t 

Figure 5b. Same as figure 5a, except here for the value of w time-averaged from to t. Early on, when radiation dominates the equation of state, 
(w) is approximately +1/3, but it tapers off and eventually reaches the value -1/3 at the present time. 



tion for w and Rh(t), so these are calculated numerically, along with R r (t). For completeness, we 
show the expansion factor a(t) as a function of cosmic time in figure 5a; the value of w, averaged 
over time from to t, is shown as a function of cosmic time in figure 5b; the gravitational radius 
R h is shown as a function of t in figure 5c; and the null geodesic itself is plotted in figure 5d. The 
blending of three different components in p produces some discernible differences in the behavior 
of R y compared to the cases we considered previously, but even here, its maximum value was a 
fraction ~0.45 of today's gravitational radius Rh(t ) at roughly 0.3t after the big bang. 

In all these cases, we have demonstrated through direct computation that null geodesies have 
never traversed the current gravitational (or Hubble) radius, regardless of what specific cosmology 
one assumes. In order for light to reach us today, its source radiating at cosmic time t e must 
have been located a proper distance R{t e ) away, corresponding to one of the trajectories shown 
in figures 1,2,3,4, or 5d. In other words, none of the sources whose light we detect today could 
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Figure 5c. Same as figure 5a, except here for the gravitational radius Rh as a function of t. 

1 . I . . . 1 . . . 1 . . . 1 . . . 1 . . . 1 

0.8 - 



a o.6 - 




0.0 0.2 0.4 0.6 0.8 1.0 



t/t 

Figure 5d. The null geodesic as a function of cosmic time t for ACDM. 

have been more distant than (in most cases, actually only a fraction ~0.3-0.4 of) our current 
gravitational horizon Rh(to)- 

In spite of this, it is sometimes claimed (see, e.g., van Oirschot, Kwan and Lewis 2010) that 
we see sources from beyond the Hubble radius because its redshift is only ~2, much smaller than 
sources such as the recombination region that produced the CMB (at a redshift of ~1 100 in the 
standard model). This misconception arises because — written as a proper distance (see Melia & 
Shevchuk 201 1) — R h is expressed as 

R h (t) = a(t)r h , (10) 
where the comoving radius may further be written in the form 

J te a(t') 

So, for example, if we crudely assume that the Universe's expansion has been driven by a matter- 
dominated p, for which a{t) = t 2/3 , we find that 

R h (t Q ) * 3ct - 3ct 2/ V /3 , (12) 
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and using the observation that Rh(t ) ~ ct (Melia & Shevchuk 201 1), with 

l +Z = ^ (13) 
a(t e ) 

(see, e.g., Weinberg 1972), one easily "finds" a redshift z h ~1.25 for the Hubble sphere (the value 
is closer to 2 when using the actual form of a(t) from ACDM). 

The apparent conflict between this inference and the curves shown in figures 1,2,3,4, and 5d, 
arises from an incorrect interpretation of R h (t) as a null geodesic, which it is not. In other words, 
the comoving radius r h in Equation (11), and therefore the proper radius R^, do not track the 
propagation of light through the Hubble flow. Only Equation (8) for R y (t) can do this because 
it includes both the effects of an expanding medium and the change in comoving radius r as 
light approaches us. Instead, Rhih) represents today's proper distance to sources that in the past 
radiated the light we see at to redshifted by an amount Zh- Therefore, to correctly determine their 
proper distance at the time the light was produced, one needs to calculate R h (t e ), not Rh(t ), where 
t e is the cosmic time corresponding to redshift z h . Clearly, z h is that special redshift for which 
Rh(t e ) = Ry{t e ), and there is therefore no conflict between the finite value of z h and the fact that 
Rh(to) is the limit to our observability today. 



5 CONCLUDING REMARKS 

Throughout this paper, we have made a conscious effort to discuss the properties of null geodesies 
in FRW cosmologies without resorting to conformal diagrams. This approach, also used to great 
effect by Ellis and Rothman (1993), makes it easier to think in terms of familiar quantities (proper 
distances and proper time) that are not always straightforward to interpret otherwise. Misconcep- 
tions often arise from the misinterpretation of coordinate-dependent effects. In their paper, Ellis 
and Rothman clearly delineated true horizons from apparent horizons, and extended the defini- 
tions, first introduced by Rindler (1956), in a clear and pedagogical manner. In this paper, we 
have paid particular attention to the gravitational horizon, also manifested as the Hubble radius, 
which is time-dependent and may or may not turn into an event horizon in the asymptotic future, 
depending on the equation of state. 

It is useful at this point to be absolutely clear about how far light has traveled in reaching 
us. It is quite evident from our results that sources whose light we see today were at a proper 
distance R(t e ) — » when t e — » 0. For the FRW cosmologies we have considered here (which 
include ACDM) the light reaching us today — including that from the recombination region asso- 
ciated with the CMB — has traveled a net proper distance of at most ~0.3-0.4c% It is therefore not 
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correct to claim that the size of the visible universe in these cosmologies is ct (or even greater in 
some interpretations). Because all causally connected sources in an expanding universe began in 
a vanishingly small volume as t —> 0, the maximum proper distance from which we receive light 
today must necessarily be less than ct , since presumably there were no pre-existing sources at a 
non-zero proper distance prior to t = with which we were in causal contact. 

It is true, however, that more and more sources become visible to us as time advances, since 
for t > to, the geodesic curves terminating in our future all rise above their current counterparts 
shown in figures 1,2,3,4, and 5d. In our future, we will see light from sources that radiated at 
proper distances greater than those shown here. Of course, Rh(t) will also continue to increase, and 
it is not difficult to convince oneself from Equations (5) and (8) that the limits of observability will 
always be R h (t), since R y (t e )/R h (t) is smaller than 1 for all t e < t. To see this, let us first find the 
emission time t e , max at which R 7 (t e ) attains its maximum value. Clearly, this happens when R y = 0, 
which means that R y (t e ^ max ) = Rh(t e , max)- But for all w > -1, Equation (5) shows that R h > 0, and 
therefore R h (t) > Rh(t e , max) for t > t e , max , which also means that R h (t) > R y {t e ^ max ) for all t > t e< max . 
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